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1. Introduction 

Experimental measurements of CP violation can be used to extract information about the CKM 
matrix. In particular, the size of indirect CP violation in the neutral kaon system, % [Q], combined 
with theoretical input, places an important constraint on the apex of the CKM unitarity triangle [^]. 
The dominant source of error in this procedure is the uncertainty in the lattice determination of the 
nonperturbative constant Bk , which parameterizes the ^° — K° matrix element. 

Promising calculations with partially quenched staggered fermions are in progress [f|, ||]. Be- 
cause staggered fermions are computationally cheaper than other standard discretizations, they 
allow simulations with the lightest dynamical quark masses currently available. Unfortunately, 
staggered fermions come with their own source of error - taste symmetry breaking. Each lattice 
fermion flavor comes in four tastes, which are degenerate in the continuum. The nonzero lattice 
spacing, a, breaks the continuum taste symmetry at 0(a 2 ), and the resulting discretization errors 
are numerically significant at present lattice spacings Thus one must use the chiral perturba- 
tion theory functional forms for staggered fermions in order to correctly perform the combined 
continuum and chiral extrapolations incorporating taste violations [f7], |h. 

2. Bk with Staggered Quarks 

The kaon B-parameter is defined as a ratio of matrix elements: 

JC K = {K°\0 K \K°) = B K J^ , (2.1) 
where &k is a weak operator and ^# vac is the result for ^k in the vacuum saturation approximation: 
6k = [s^Yn (1 + Ys)d a } [sEYn ( 1 + Ys )4] (2-2) 
^vac = ^nsaY^ + Y5)da\\0)(0\[sbY f i^ + Y5)db}\K°). (2.3) 

Note that there are separate summations over color indices a and b. Thus the matrix element jft^ 
receives contributions from two different quark-level contractions, one of which produces a single 
color loop and the other which has two color loops. 

In order to calculate Bk with staggered fermions we must introduce the taste degree of free- 
dom, both in the operators an d in the states (K° and K°). We choose the external staggered 
kaons to be taste P, by which we mean that the lattice meson operator contains the pseudoscalar 
taste matrix £5. Since this is the lattice Goldstone taste, its correlation functions satisfy U(\)a 
Ward identities, so S^PT expressions for its mass, decay constant, and other physical quantities 
are simpler than those for other tastes of PGBs. In addition, it is a local kaon on the lattice, and 
therefore relatively simple to implement numerically. Because the external kaons are taste P, the 
weak operator should also be taste P: 

@T e = [SafoCl + 7s) ® $s)da] [*(y M (l + Ys) ® (2-4) 

However, when this operator is Fierz-transformed, it mixes with other tastes. To remedy this, we 
introduce two types of valence quarks, 1 and 2, into the Bk operator: 

fpggered = 2 { [^(y^ ( 1 + 7s ) ® ^ 5 )jl fl ] ^(/ M (l + y 5 ) ® $ 5 )d2 b ] 

+ + 0$) \^b{Y^ + Y5)®^)d2 a ]\ (2.5) 
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and take the matrix element between two types of kaons: 

^staggered = ^ | ^aggered ^ ^ (2 g) 

where K® is an sldl meson and K$ is an s2d2 meson. The extra valence quarks require explicit 
inclusion of both color contractions in 0^ aggered > while the overall factor of 2 ensures that ^k and 
^staggered ^ aye same total number of contractions. When ff^" ggered [ s Fierz-transformed, it 
has a new flavor structure so it does not contribute to ^ K taggered . Thus ^ K taggered cannot receive 
contributions from incorrect tastes. 

Current calculations of Bk with staggered quarks use the partially quenched approximation. 
The quark content the corresponding PQ theory is quite large. We have already seen that staggering 
requires two sets of d and s valence quarks. Partial quenching adds two corresponding sets of ghost 
quarks as well as three sea quarks, resulting in eleven total quark flavors, each of which comes in 
four tastes. In order to make this completely clear, we show the explicit form of the quark mass 
matrix: 

M = diag{m u I, m^I, m s I, m x I, m y I, m x I, m y I, m x I, m y I, m x I, m y I}, (2.7) 

sea valence 1 valence 2 ghost 1 ghost 2 

where / is the 4x4 identity matrix. 

3. Generalized S^PT for B K 

The expression for Bk in S^PT must describe Bk at a ^ if it is to be used for continuum 
and chiral extrapolations of lattice data; we therefore discuss how Bk is actually calculated on the 
lattice. Additional details can be found in Ref. [jlj]. The Bk matrix element, ^ K taggered , receives a 
contribution from the lattice version of £^ ag8er , as well as from other lattice operators that are in 
the same representation of the symmetry group that maps a hypercube onto itself ^ : 

^aggered,cont = ^taggeredjat _« , fc p , U_ ^ 

K K An An 

+a 2 [all taste ops.] +a ,2 [all taste ops.] + . . . , (3.1) 

where a is the strong coupling constant. The 1-loop perturbative matching coefficients between 
^staggere ^ ^ e contmuum anc j four-fermion lattice operators are known, and are numbers of order 



unity times ajAn [10]. However, the 2-loop matching coefficients have not been determined, so, 
in order to remain conservative, we consider them to be of order unity times a 2 without any factors 
of An. Current numerical staggered calculations are in fact of the following matrix element: 

{K° p \ ^_ /oop (taste P) \K* P ) = ^f lat , (3.2) 

where the subscript "1 — loop" and the argument "taste P" indicate that one includes all staggered 
lattice operators with taste P that mix with the latticized Bk operator at GiajAn), i.e. those in 
the second term on the RHS of Eq. ( |3.1[ ), using the appropriate matching coefficients. However, 
one neglects wrong-taste and higher-order perturbative mixing (terms three and four), as well as 



all operators which arise through discretization effects (term five), in Eq. (3.1). Although the 
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expression in Eq. ( ^2| ) differs from the continuum matrix element, it reduces to the desired quantity 
in the continuum limit. Generically, 

(X 

= .^cont + -7-Jt' + a 2 J(" + a 2 J("' + ..., (3.3) 
An 

where ^# con t is the desired continuum result. The matrix element comes from neglecting taste- 
violating 1-loop operator mixing, while ^M" comes from neglecting 2-loop operator mixing. Both 
taste-breaking and taste-conserving discretization errors generate M'" . Our goal is therefore to 
calculate the lattice matrix element, ^#i a t, to NLO in S^PT. One can fit the resulting expression 
to lattice data, determine the various errors, and remove them, thereby extracting the continuum 
matrix element, ^# con t- 

The appropriate power-counting scheme for calculating in S^PT must incorporate a 2 , 
a /An, and a 2 . Continuum #PT is a low-energy expansion in both the pseudo-Goldstone boson 
(PGB) momentum and the quark masses; it assumes that p 2 PGB /A 2 x ~ m 2 p GB /K 2 x , where mp GB oc m q 
and A x « Anf % is the ^PT scale. However, the numerical values of m q , a and a all depend on 
the particular parameters of a given lattice simulation. Current PQ staggered lattice simulations 
[|]] use a range of PGB masses from tnp GB /Al = 0.04 — 0.2, so our S^PT expression must apply 
throughout this range. Generic discretization errors are of the size g! 2 Aq CD , which is approximately 
0.04 for 1/a ~ 2GeV and Aqcd ~ 400 MeV, so they are comparable to the minimal mp^/A 2 ^ 
and should be included at the same order. Taste-breaking discretization errors, on the other hand, 
are caused by exchange of gluons with momentum n/a, and therefore receive an additional factor 
of Oy(n/a), so their size must be considered separately. At the lightest quark masses, the lattice 
Goldstone pion mass is comparable to the mass splittings among the other PGB tastes: vn 2 PGB l h 2 % ~ 
a 2 Oy(q* = 7i/a)A 2 ~ 0.04, where A is a QCD scale which turns out to be around 1200MeV. 
Thus they are not suppressed relative to pure Gia 2 ) discretization effects by the additional powers 
of a, as naive power-counting would suggest, because of the large scale A associated with the 
taste-breaking process at the quark level Standard S^PT only includes discretization effects 
- we now consider additional errors from perturbative operator matching, which depend upon 
0Cy (<?*). Generically, the choice of q* is process-dependent, and the value of ay{q*) ranges from 
~ 0.3 — 0.55 for q* = n/a — I /a at a = .125fm [Q. Thus both a /An and a 2 must be included at 
lowest-order in our power counting. 

In light of this discussion, we adopt the following extended S^PT power-counting scheme: 

p 1 ~ m ~ a 2 ~ a 2 a ~ a/An ~ a 2 , (3.4) 

where a 2 a = oty (%/a)a 2 . We account for the fact that a 2 a terms in the action are enhanced nu- 
merically by including them at the same order as simple discretization effects. In fact, while it 
may seem ad hoc, p 2 ~ m ~ a a is the standard S#PT power counting scheme. It is simply not 
traditionally written as such because standard S^PT calculations have only included &{a 2 a ) taste- 
breaking discretization errors from the action, and have therefore not needed to contrast them with 
pure &{a 2 ) discretization effects. We also use conservative power-counting for the perturbative 
errors by assuming that 2-loop contributions are not significantly smaller than those from 1-loop 
diagrams. We emphasize that our scheme is phenomenologically based on the particular parameter 
values of current staggered simulations - simulations using significantly lighter quark masses or 
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smaller lattice spacings would require a different scheme and result in a different S^PT expression 
for B K . 

In S^PT, ^#i at is simply the matrix element of a sum of operators, 

^lat = (^ P |^4l*2/>>> < 3 " 5 ) 

with undetermined coefficients. In order to calculate ^#i at at NLO we must include operators in the 
chiral effective theory of the following order: 



-), 0(m), 0(a z ), #{a z a ), 0(a/4it), 6(a 2 ). (3.6) 
We determine all of the linearly-independent operators of this order using the graded group the- 



ory method of Ref. pi3[ ] along with the taste spurion method of Refs. [|I4|, m. Directly mapping 
Q^taggeied ontQ c hiral operators leads to operators of of &(p 2 ), truncation of perturbative matching 
factors leads to operators of 6{ajA%),G(a 2 \ generic discretization effects lead to operators of 
&(a 2 ), and insertions of taste-breaking operators from the staggered action with o K agsered lead to 
operators of &{a 2 a ). In total, we find thirteen chiral operators; they are given in Eq. (27) of Ref. |jl|]. 
We note, however, that many of these chiral operators arise in more than one way and correspond to 
more than one quark level operator, so they actually have more than one undetermined coefficient. 

4. Results and Conclusions 

We have calculated Bk to NLO in S^PT for quenched, PQ, and full QCD; the results are in 
Ref. [[lj]. The simplest expression that still shows the general structure of Bk in S^PT is the partially 
quenched expression for for three sea quarks and degenerate valence quarks (m x = m y ): 



The functions £ and i are chiral logarithms, and the index B' runs over the PGB tastes I,P,V,A, 
and T. The fact that all PGB tastes enter through loops is a generic trait of S^PT. The 1-loop 
contribution in Eq. ( fOl ) proportional to Bo and the analytic terms proportional to the coefficients B 
and C reproduce the continuum #PT result when a — > [|5). The 1-loop contributions proportional 
to the coefficients 'ta lB and are new, and come perturbative matching and discretization eiTors. 
Because they have a different functional form from the continuum 1-loop expression, they can, in 
principle, be determined and removed separately at each lattice spacing. There remains, however, 
a multiplicative correction to Bo from the discretization and perturbative matching errors, A, which 
can only be removed using a fit to multiple lattice spacings. 

Because of the large number of fit parameters and complexity of the expression, we outline 
a fitting strategy that takes advantage of other matrix elements and maximizes information that 
can be extracted from a single lattice spacing in Ref. Even with the implementation of these 
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suggestions, however, fitting the staggered Bx data will be highly nontrivial. One would clearly like 
to somehow decrease the number of operators, or at least undetermined coefficients, that contribute 
to Bk at NLO. One way to do this is by using improved links. This drastically reduces the 1- 
loop perturbative mixing between &\-i 00 p (taste P) and certain wrong-taste operators. The biggest 
source of potential improvement, however, is in better perturbative matching. A good first step 
would therefore be to match to all tastes of lattice operators at 1-loop, since the requisite matching 
coefficients are known. While it has been thought that "wrong taste" operator contributions to 
Bk would be less important than those from taste P operators, our power-counting and operator 
enumeration show that this is not the case in general. Fully nonperturbative matching, although 
extremely difficult, would completely eliminate some of the thirteen operators. 

Bk is an important parameter for constraining the phase of the CKM matrix, and consequently 
new physics. Calculations of Bk with staggered fermions are promising, but they need the appro- 
priate S^PT form for correct continuum and chiral extrapolations. We have calculated Bk to NLO 
in S^PT for quenched, PQ, and full QCD using an extended power-counting scheme. The result is 
more complicated than in previous S^PT cases. While in expressions for quantities such as fx and 
fo, the primary effect of staggering is additive corrections to PGB masses inside loops, the expres- 
sion for Bx contains entirely new contributions with different functional forms than the continuum 
piece. Use of our expression for Bx at NLO in S#PT, in combination with sufficient lattice data, 
should allow a precise determination of Bx with staggered quarks. 
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